Dynamic dissipative cooling of a mechanical oscillator in strong-coupling 
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Cooling of mesoscopic mechanical resonators represents a primary concern in cavity optomechan- 
ics. Here in the strong optomechanical coupling regime, we propose to dynamically control the 
cavity dissipation, which is able to significantly accelerate the cooling process while strongly sup- 
pressing the heating noise. Furthermore, the dynamic control is capable of overcoming quantum 
backaction and reducing the cooling limit by several orders of magnitude. The dynamic dissipation 
control provides new insights for tailoring the optomechanical interaction and offers the prospect of 
exploring macroscopic quantum physics. 
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One of the ultimate goals in quantum physics is to over- 
come the thermal noise, so that quantum effects can be 
observed experimentally. A prominent example is cavity 
optomechanics [l|, Q , which enables not only the funda- 
mentally test of quantum theory and the exploration of 
quantum-classical boundary, but also important applica- 
tions in quantum information processing and precision 
metrology. For these applications, the first crucial step 
is to prepare the mechanical resonator into the quantum 
regime So far, numerous experiments have focused 

on backaction cooling in the weak optomechan- 

ical coupling regime, holding potential for ground-state 
preparation of mechanical resonators in the resolved side- 
band condition [13Hl5l |. along with backaction evading 
quantum non-demolition measurements Further 
step lies in the strong coupling, essential for coherent 
quantum optomechanical manipulations and 
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electromechanical interactions [25|, [26|. However, till 
date, strongly-coupled optomechanical cooling has pre- 
dicted only limited improvement over weak coupling due 
to the saturation effect of the steady-state cooling rate 
22, 27, 28]. Although strong coupling allows state swap- 
ping [g, |21| , it cools the mechanical resonator only at a 
single instant in the Rabi oscillation cycle. Thus it is ur- 
gent to overcome these limitations for cooling and manip- 
ulating mesoscopic mechanical systems in the quantum 
regime. 

For this purpose, in this Letter we show the dynamic 
tailoring of the cooling and heating processes by exploit- 
ing the cavity dissipation, overcoming the saturation of 
the steady-state cooling rate. This greatly accelerates the 
cooling process and thereby strongly suppresses the ther- 
mal noise. Moreover, heating induced by swapping and 
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FIG. 1: (color) (a) Sketch of a typical optomechanical system, 
(b) Level diagram of the linearized Hamiltonian fTJ. \n,m) 
denotes the state of n photons and m phonons in the displaced 
frame. The solid (dashed) curves with arrows correspond to 
the cooling (heating) processes. 



interaction quantum backaction are largely suppressed 
by periodic modulation of the cavity dissipation, which 
breaks the fundamental limitation of backaction cooling. 

We consider a generic optomechanical system in which 
an optical cavity driven by a laser is coupled to a me- 
chanical resonance mode, as illustrated in Fig. [TJa). 
In the rotating frame at the driven laser frequency ui, 
the system Hamiltonian reads H = — (w — u c )aJa + 
w m tfb + ga^a(b + &t) + (Oa* +Q,*a) where a (b) 

represents the annihilation operator for the optical (me- 
chanical) mode with uj c (cj m ) being the corresponding an- 
gular resonance frequency; g denotes the single-photon 
optomechanical coupling rate; fl represents the driving 
strength. For strong driving, the Hamiltonian can be 
linearized, with a = a\ + a, b = b± + j3. Here a\ and b\ 
describe the fluctuations around the mean values a = (a) 
and j3 = (b), respectively. Neglecting the nonlinear 
terms, this yields the Hamiltonian 

H L = -A'a t 1 a 1 + ujt^^ + (Ga\ + G*a x ){b l + b{), (1) 

where A' = u — lj c + 2|G7| 2 /cj m is the optomechanical- 
coupling modified detuning, and G — ag describes 
the linear coupling strength. Taking the dissipa- 
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tions into consideration, the system is governed by the 
quantum master equation p = i[p, Hi] + re2?[ai]p + 
7 (n t h + l)f [h]p + in th T>[b\]p, where X>[o]p = <5po t - 
(o^op + p6 t 6)/2 denotes the Liouvillian in Lindblad form 
for operator 6; re = lo c /Q c (7 = w m /Q m ) represents 
the dissipation rate of the optical cavity (mechanical) 
mode; nth = l/(e fiWm/,fc B T — 1) corresponds to the ther- 
mal phonon number at the environmental temperature 
T. 

Figure QJb) displays the level diagram of Hi and the 
coupling routes among states \n, to) with n (to) being 
the photon (phonon) number in the displaced frame. We 
note that there are three kinds of heating processes de- 
noted by the dashed curves in Fig. Hlb), corresponding to 
swap heating (_B), quantum backaction heating (D) and 
thermal heating (F). Suppressing thermal heating is the 
ultimate goal while swap heating and quantum backac- 
tion heating are the accompanying effect when radiation 
pressure is utilized to cool the mechanical motion. Swap 
heating emerges when the system is in the strong cou- 
pling regime which enables reversible energy exchange 
between photons and phonons. Meanwhile, quantum 
backaction heating can pose a fundamental limit for back- 
action cooling. The solid curves (A, C and E) illus- 
trate cooling processes associated with energy swapping, 
counter-rotating- wave interaction and cavity dissipation, 
which one seeks to enhance while suppressing heating for 
efficient mechanical motion cooling. 

We focus on the resolved sideband regime re < w m and 



we set A' = — l 



in which the beam splitter interaction 



a\b\ + a\b\ is on resonance. In this case the dynamical 
stability condition from the Routh-Hurwitz criterion [3(il ] 
requires 2 \G\ < oj m . To realize cooling, the cooperativity 
C = 4|G| 2 /(7k) 3> 1 should also be satisfied. Starting 
from the master equation, we obtain a set of differential 
equations for the mean values of the second-order mo- 



ments N a = (a|a 1 ), Nf, — (6{6 1 ), (al&i), (ai&i), (a 2 ) and 
(b\) (see Supplementary Material [3l| ) . In the steady 



state we obtain the phonon occupancy 22, 2 
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where the first term is the classical cooling limit and 
the second term originates from the quantum backac- 
tion, consisting of both dissipation quantum backaction 
related to the cavity dissipation (with the associated 
fluctuation-dissipation theorem) and interaction quan- 
tum backaction associated with the optomechanical in- 
teraction (see Supplementary Material 3l| for full de- 
scription). In the weak coupling regime, Eq. (|2|) reduces 

to ~ 7 «th/(r + 1) +^1(1^1) with r = 4|G| 2 /«, 



which agrees with Refs. [13|, |14(, and with re 2 / (16^) 
the dissipation quantum backaction from fluctuation- 
dissipation. In the strong coupling regime, we obtain 
|Gf/[2(u£ - 4|G| 2 )]. In this case 
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FIG. 2: (color) (a) Time evolution of mean phonon number 
N b for G/u m = 0.005, 0.01, 0.02 and 0.1 (numerical results), 
(b) Nb for G/u)m = 0.005 and 0.01 with a wider time inter- 
val. The shadowed region shows the same time interval with 
(a), (c) Modulation scheme of the cavity dissipation rate n(t) 
for fast cooling to the steady-state limit and (d) the time 
evolution of mean phonon number Nt with (red solid curve) 
and without (blue dashed curve) modulation for G/w m = 0.2. 
Other parameters: n t h = 10 3 , 7/u; m = 10~ 5 , ft/oj m = 0.05. 
The dotted horizontal lines correspond to the steady-state 
cooling limits, given by Eq. @. 



the classical limit is restricted by the cavity dissipation 
rate re, while the interaction quantum backaction limit 
suffers from high coupling rate |G|. 

To study the cooling dynamics beyond the steady 
state, we solve the differential equations to obtain the 
time evolution of the mean phonon number N b . For 
weak coupling, we have iV h wk ~ 7i th (7 + re~ rt )/(7 + T) + 
[re 2 /(16w 2 n )](l — e _rt ), which shows that the mean 
phonon number decays exponentially with the cooling 
rate T. This cooling rate is limited by the coupling 
strength, since in the cooling route A — > E the energy 
flow from the mechanical mode to the optical mode (pro- 
cess A) is slower than the cavity dissipation (process E) . 

In the strong coupling regime, we obtain the time evo- 
lution of the mean phonon number described by (see Sup- 
plementary Material |3lj ]) 
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where lj± = lo^ ± 2|G|w m are the normal eigenmode 
frequencies. The phonon occupancy exhibits oscillation 
under an exponentially-decaying envelope and can be di- 
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vided into two distinguished parts and N^, where 
the hrst part originates from energy exchange between 
optical and mechanical modes, and the second part is 
induced by quantum backaction. reveals Rabi oscil- 
lation with frequency ~ 2\G\, whereas the envelopes have 
the same exponential decay rate V = (k + j)/2 regard- 
less of the coupling strength \G\. This is because, in the 
strong coupling regime, the cooling route A — > E is sub- 
jected to the cavity dissipation (process E), which has 
slower rate than the energy exchange between phonons 
and photons (process A). This saturation prevents a 
higher cooling speed for stronger coupling. In Figs. HJa) 
and (b) we plot the numerical results based on the master 
equation for various G. It shows that for weak coupling 
the cooling rate increases rapidly as the coupling strength 
increases, whereas for strong coupling the envelope decay 
no longer increases, instead the oscillation frequency be- 
comes larger. 

Fast cooling to the steady-state limit. — To speed up 
the cooling process in the strong coupling regime, here 
we take advantage of high cavity dissipation to dynami- 
cally strengthen the cooling process E. The internal cav- 
ity dissipation is abruptly increased each time when the 
Rabi oscillation reaches a minimum-phonon state, such 
as through RF-synchronized carrier injection to the opti- 
cal cavity [32| • At this time the system has transited from 
state \n,m) to state \n + l,m — 1). Once a strong dissipa- 
tion pulse is applied to the cavity so that the process E 
dominates, the system will irreversibly transit from state 
| n + 1,to — 1} to state 1 7i, 777 — 1). The dissipation pulse 
has essentially behaves as a switch to halt the reversible 
Rabi oscillation, resulting in the suppression of the swap 
heating. To verify this dissipative cooling, in Figs. [2] (c) 
and (d) we plot the modulation scheme and the corre- 
sponding time evolution of mean phonon number Nb for 
rt/wm = 0.05 and G/w m = 0.2. At the end of the first half 
Rabi oscillation cycle, t ~ ir/(2\G\), a dissipation pulse of 
pulscwidth 0.0l7r/(2|G|) is applied. Detailed tradeoffs of 
the dissipation quantum backaction and the interaction 
quantum backaction for varying _pulsewidths are shown 
in the Supplementary Material [3l|. After incidence of 
the dissipation pulse, the phonon number reaches and re- 
mains near the steady-state limit. For short time scales, 
the remaining small-amplitude oscillations mainly origi- 
nate from counter-rotating-wave interactions. Note that 
without modulation (blue dashed curve), the steady-state 
cooling limit is reached only after t ~ 400/u; m ; while with 
the modulation (red solid curve), it only takes t ~ 8/w m 
to cool below the same limit. 

Breaking the fundamental limit of backaction cooling. 
— By periodically modulating the cavity dissipation so as 
to continuously suppress the swap heating, the phonon 
occupancy can be kept below the steady-state cooling 
limit. Actually, each time after the dissipation pulse is 
applied, the photon number quickly drops to the vac- 
uum state, which equivalently re-initializes the system. 
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FIG. 3: (color) Modulation scheme of n(t)/n(0) (a) and 
the corresponding time evolution of Nb (b) for G/ui m = 
0.1, K(0)/w m = 0.01 (red solid curve) and 0.02 (blue dashed 
curve). In (b), the two dotted horizontal lines (from top to 
bottom) denoting the respective steady-state cooling limits 
depending on the cavity decay k(0), given by Eq. the 
dash-dotted line denotes the instantaneous-state cooling limit 
independent of k(0), given by Eq. ©; the "ON" and "OFF" 
regions corresponds that the modulation is turned on and off, 
respectively; the vertical coordinate range from 10 to 10 3 is 
not shown. Other parameters: n t h = 10 3 , j/u m = 10~ 5 . 



By periodic pulse application, the system will periodi- 
cally re-initializes, which keeps the phonon occupancy in 
an instantaneous-state cooling limit as verified in Fig. 
[3] The instantaneous-state cooling limit is given by (see 
Supplementary Material (3lj ) 
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Here the first term comes from for t ~ 7r/(2|G|), 
which shows a 7tk/(4|G|) times reduction of classi- 
cal steady-state cooling limit. The second term of ~ 
7r 2 |G| 4 /uif n , obtained from when t ~ n/uj m , reveals 
that the second order term of |G|/o> m in quantum back- 
action has been removed in our approach, leaving only 
the higher-order terms. Note that the cooling limit (|4| is 
the sum of the individual minimum of iVJ ^ and % in 
their first oscillation cycle. Notably, in Fig. [3] we demon- 
strate that the modulation is switchable. If we turn on 
the modulation ("ON" region), the system will reach the 
instantaneous-state cooling limit; if we turn off the mod- 
ulation ("OFF" region), the system transits back to the 
steady-state cooling limit. 

In particular, from Eq. Q, the interaction quantum 
backaction heating term forms a carrier-envelope 

type evolution, where the carrier oscillation represents 
the counter-rotating-wave interaction and the envelope 
oscillation is a result of coherent energy exchange due 
to strong coupling. The minimum of A|*| is depen- 
dent on the carrier-envelope frequency matching. If 
(lu + + uj-)/(uj + — oj-) — k (k — 3,5...), yielding 

reaches a mini- 



|G|/w m = k/(k z + 1) = 0.3, 5/26.. 
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FIG. 4: (color) Different pulsewidth (a) and the correspond- 
ing time evolution (b) of jV&. The horizontal lines indicate 
the three cooling limits given by Eqs. (J5J), (© and ^ from 
top to bottom. The inset shows the effective cooling rate as 
a function of time. Parameters : nth = 300, G/w m = 0.3, 
K./com = 0.003 and y/u m = 10" 5 . 
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for t ~ 7r/ (2|G|). Here we obtain the 
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optimized instantaneous-state cooling limit as (|3l|) 
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which reduces both the classical and quantum steady- 
state cooling limits by a factor of 7tk/(4|G|). Remarkably, 
this reduction is significant when the system is in the 
deep strong coupling regime. Besides, the leading order 
of the interaction quantum backaction heating scales as 
k\G\/ui^ which can be a few orders of magnitude lower 
than the steady-state case, representing large suppression 
of quantum backaction. 

To verify suppression of the interaction quantum back- 
action heating, in Fig. |4] we plot the cooling dynam- 
ics with dissipation modulation for G/u) m = 0.3 and 
n/uj m — 0.003. The single modulation pulse brings down 
the phonon occupation to the optimized instantaneous- 
state cooling limit described in Eq. ([5]), with the time- 
dependent effective cooling rate T e g = (dNb/dt)/N~b 
shown in the inset. With short-pulse modulation (blue 
dashed curve), the remaining oscillation, mainly induced 
by the counter-rotating-wave interaction, has a quasi- 
period of 7r/(2|G|) due to frequency matching. This 
small-amplitude fluctuations around the instantaneous- 
state cooling limit might affect future quantum proto- 
cols, but in the sense of time-averaged through timescales 
larger than 7r/(2|G|), the cooling limit can be viewed sta- 
ble. By using long-pulse modulation, the quasi-periodic 
fluctuations can be suppressed (red solid curve). This is 
because the large dissipation suppresses the interaction 
quantum backaction. The cost is that the dissipation 
quantum backaction takes effect and gradually increases 
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FIG. 5: (color) Cooling limits given by Eqs. ([2]) (black 
dotted curves), (j4]) (blue dashed curves) and ([5j (red solid 
curves) versus G/w m for n t h = 10 3 (thin curves), 3 x 10 2 
(thick curves). Other parameters are n/uj m = 0.003 and 
7 /w m = 10" 5 . 



the phonon number. This trade-off can be balanced by 
optimizing the pulsewidth as shown in the Supplemen- 
tary Material [3 If . 

Figure [5] plots the cooling limits as functions of G/u) m , 
which reveals that instantaneous-state cooling limits are 
much lower than steady-state cooling limits. For small 
coupling rates, we observed that interaction quantum 
backaction is insignificant and suppression of swap heat- 
ing is the main origin of cooling limit reduction. For large 
coupling rates, suppressing interaction quantum backac- 
tion is crucial for obtaining lower limits. Typically, the 
cooling limits can be reduced by a few orders of magni- 
tude. For example, when G/w m = 0.3 and n/uj m — 0.003, 
we obtain N st & = 3.4, while = 0.03, corresponding 

to more than 100 times of phonon number suppression. 

Experimentally, the dynamic control of cavity dis- 
sipation can be realized, for example, by modulating 
free-carrier plasma density [32T434I or using a light ab- 
sorber /scatterer [3a |. Note that we assume G is kept un- 
changed when the dissipation pulses are applied, which 
corresponds to the invariablencss of the intracavity field 
a. This can be fulfilled by simultaneously changing the 
driving f2(t), so that equation [iA' — K(t)/2]a — iCl(t) = 
is satisfied all the time ([3l|; Section VI). Here modulated 
square-shaped dissipation pulses are used though further 
simulations show that the results are irrespective of the 
pulse shape, as long as they are executed quickly with 
strong enough peak value at the desired time. This is 
because the pulse dissipation mainly relies on the pulse 
area. 

In summary, we examined cooling of mesoscopic me- 
chanical resonators in the strong coupling regime and 
propose dynamic dissipative schemes which possess large 
cooling rates, low cooling limits, and long-time stability. 
By making use of the cavity dissipation, swap heating can 
be strongly avoided and the interaction quantum backac- 
tion largely suppressed, with great advantages over the 
current conventional cooling approaches. For example, 
a single dissipation pulse enables more than 50 times 
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higher cooling rate; with periodic modulation of cavity 
dissipation, the cooling limit can be reduced by more 
than two orders of magnitude. Different from the cool- 
ing schemes with modulated coupling [3(il - l40j |. we take 
advantage of large cavity dissipation, usually regarded 
as a noise source. Together with recent proposals of 
other dissipative effects such as two-level ensembles 41 1 
or photothermal effect 42 , 43| , we demonstrate that cav- 
ity dissipation (even in the presence of the considered 
dissipation quantum backaction) can be viewed as a re- 
source. Compared with the dissipative coupling (44 - 46 1 . 
this active dissipation control does not require the cou- 
pling between the cavity dissipation and the mechanical 
resonator. The dynamic dissipative cooling provides a 
new way for exploring the quantum regime of mechani- 
cal devices, ranging from mechanical ground state prepa- 
ration, to generation of mesoscopic quantum states, and 
quantum- limited measurements. 
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